Geographical networks evolving with an optimal policy 
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In this article, we propose a growing network model based on an optimal policy involving both 
topological and geographical measures. In this model, at each time step, a new node, having 
randomly assigned coordinates in a 1 x 1 square, is added and connected to a previously existing 
node i, which minimizes the quantity r 2 /kf, where r, is the geographical distance, ki the degree, and 
a a free parameter. The degree distribution obeys a power-law form when a = 1, and an exponential 
form when a = 0. When a is in the interval (0, 1), the network exhibits a stretched exponential 
distribution. We prove that the average topological distance increases in a logarithmic scale of 
the network size, indicating the existence of the small-world property. Furthermore, we obtain the 
geographical edge-length distribution, the total geographical length of all edges, and the average 
geographical distance of the whole network. Interestingly, we found that the total edge-length will 
sharply increase when a exceeds the critical value a c = 1, and the average geographical distance 
has an upper bound independent of the network size. All the results are obtained analytically with 
some reasonable approximations, which are well verified by simulations. 

PACS numbers: 89.75.Hc, 87.23.Ge, 05.40.-a, 05.90.+m 



INTRODUCTION 



Since the seminal works on the small-world phe- 
nomenon by Watts and Strogatz |1] and the scale-free 
property by Barabasi and Albert [2| , the studies of com- 
plex networks have attracted a lot of interests within 
the physical community [H, 0, H, @- Most of the pre- 
vious works focus on the topological properties (i.e. non- 
geographical properties) of the networks. In this sense, 
every edge is of length 1, and the topological distance be- 
tween two nodes is simply defined as the number of edges 
along the shortest path connecting them. To ignore the 
geographical effects is reasonable for some networked sys- 
tems (e.g. food webs Q, citation networks Q, metabolic 
networks Q), where the Euclidean coordinates of nodes 
and the lengths of edges have no physical meanings. Yet, 
many real-life networkSjSuch as transportation networks 
M tlil, the Internet [jj, [H[, and power grids 0, [H[, 
have well-defined node-positions and edge-lengths. In ad- 
dition to the topologically preferential attachment intro- 
duced by Barabasi and Albert [2| , some recent works have 
demonstrated that the spatially preferential attachment 
mechanism also plays a major role in determining the 
network evolution [l(| [l?], E3] • 

Very recently, some authors have investigated the spa- 
tial structures of the so-called optimal networks fl9l . |2C| . 
An optimal network has a given size and an optimal link- 
ing pattern, and is obtained by a certain global opti- 
mization algorithm (e.g. simulated annealing) with an 
objective function involving both geographical and topo- 



logical measures. Their works provide some guidelines 
in network design. However, the majority of real net- 
works are not fixed, but grow continuously. Therefore, 
to study growing networks with an optimal policy is not 
only of theoretical interest, but also of practical signif- 
icance. In this paper, we propose a growing network 
model, in which, at each time step, one node is added 
and connected to some existing nodes according to an op- 
timal policy. The degree distribution, edge-length distri- 
bution, and topological as well as geographical distances 
are analytically obtained subject to some reasonable ap- 
proximations, which are well verified by simulations. 



II. MODEL 

Consider a square of size lxl with open boundary con- 
dition, that is, a open set (0, 1) x (0, 1) in Euclidean space 
M 2 , where "x" signifies the Cartesian product. This 
model starts with m fully connected nodes inside the 
square, all with randomly assigned coordinates. Since 
there exists m nodes initially, the discrete time steps in 
the evolution are counted ast = m+l,m + 2, Then, 
at the tth time step (t > m), a new node with randomly 
assigned coordinates is added to the network. Rank each 
previously existing node according to the following mea- 
sure: 



\Tl-Tl\ 2 /kf{t), i = 1,2,-. -,t-l, 



(1) 
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and the node having the smallest ui is arranged on the 
top. Here, each node is labelled by its entering time, ~r} 
represents the position of the tth node, ki(t) is the degree 
of the ith node at time t, and a > is a free parameter. 
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FIG. 1: Some network examples for different values of a: (a) 
a = 0, (b) a = 0.5, (c) a = 1, (d) a = 2. All the four 
networks are of size N = 300 and m = 1. 



The newly added node will connect to to existing nodes 
that have the smallest u> (i.e. on the top of the queue). 
All the simulations and analyses shown in this paper are 
restricted to the specific case m = 1, since the analytical 
approach is only valid for the tree structure with to = 1 . 
However, we have checked that all the results will not 
change qualitatively if to is not too large compared with 
the network size. 

In real geographical networks, short edges are always 
dominant since constructing long edges will cost more 
[2lj . On the other hand, connecting to the high-degree 
nodes will make the average topological distance from the 
new node to all the previous nodes shorter. These two 
ingredients are described by the numerator and denomi- 
nator of Eq. (1), respectively. In addition, the weights of 
these two ingredients are usually not equal. For example, 
the airline travellers worry more about the number of se- 
quential connections pj| . the railway travellers and car 
drivers consider more about geographical distances [l(|, 
and the bus riders often simultaneously think of both 
factors [22j |. In the present model, if a = 0, only the ge- 
ographical ingredient is taken into account. At another 
extreme, if a — > oo, the geographical effect vanishes. 

Fig. 1 shows some examples for different values of 
a. When only the geographical ingredient is considered 
(a = 0), most edges are very short and the degree distri- 
bution is very narrow. In the case of a = 0.5, the average 
geographical length of edges becomes longer, and the de- 
gree distribution becomes broader. When a = 1, the 
scale-free structure emerges and a few hub nodes govern 
the whole network evolution. As a becomes very large, 
the network becomes star-like. 
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FIG. 2: Numerical simulation about the relation between Si 
and ki. The solid line represents the assumption Si tx kf . In 
this simulation, we first generate a network of size N = 10000 
by using the present optimal growing policy. To detect Si, 
some test nodes with randomly assigned coordinates are used, 
and for each test node, following Eq. (1), the existing nodes 
with minimal u is awarded one point. The test node will be 
removed from the network after testing. The area Si for the 
ith node is approximately estimated as the ratio of the score 
node i has eventually got after all the testings to the total 
number of test nodes. The simulation shown here is obtained 
from 30000 test nodes, and the parameter a = 0.9 is fixed. 



III. DEGREE DISTRIBUTION 

At the tth time step, there are t — 1 pre-existing 
nodes. The square can be divided into t — 1 regions 
such that if a new node is fallen inside the zth region 
S i} i € {1, 2, • • • ,t - 1}, the quantity |r* - n\ 2 /kf{t) 
is minimized, thus the new node will attach an edge to 
the «th node. Since the coordinate of the new node is 
randomly selected inside the square, the probability of 
connecting with the ith node is equal to the area of Si. 
If the positions of nodes are uniformly distributed, sta- 
tistically, the area of S% is approximately proportional to 
kf{t) with a time factor h(t) as Si ~ h{t)kf{t). Fig. 2 
shows the typical simulation results, which strongly sup- 
port the valid of this assumption. Accordingly, by using 
the mean-field theory 0], an analytic solution of degree 
distribution can be obtained. However, when a > 1, 
most edges are connected to one single node (see Fig. 
Id), so analytic solution is unavailable. Here, we only 
consider the case of 1 > a > 0. Assume 



At, (2) 

where A is a constant that can be determined self- 
consistcntly. Using the continuum approximation in time 
variable t, the evolving of node z's degree reads 

= = ~Ar 



dt 



(3) 
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with the initial condition kAi) = 1. The solution is 



h(t) = F(t/i), (4) 

where 

/ 1-Q \V(l-a) 

F(x)=(l+—Hx)) . (5) 

Accordingly, the degree distribution can be obtained as 

P(k) = \ £ S(k t (t) -k)= A ^_ a _ i} . (6) 
1 i=i k a e — 

The constant A = duF a (\/u) is determined by the 
condition J kP(k)dk = (k) = 2, where (k) = 2m = 2 sig- 
nifies the average degree. The above solution is similar 
to the one obtained by using the a ppr oach of rate equa- 
tion proposed by Krapivsky et al. [23|. In addition, one 
should note that if a = 1, the mean-field theory yields a 
solution P(k) ~ fc~ 3 , which is comparable to the exact 
analytic solution P(k) = A/(k(k + l)(fc + 2)). 

Clearly, the degree distribution obeys a power-law form 
at a — 1, and an exponential form at a — 0. When a is 
in the interval (0,1), the networks display the so-called 
stretched exponential distribution [24|: For small a, the 
distribution is close to an exponential one, while for large 
a, it is close to a power law. This result is in accordance 
with the situation of transportation networks. If only the 
geographical ingredient is taken into account (e.g. road 
networks [l!|), then the degree distribution is very nar- 
row. On the contrary, if the topological ingredient plays 
a major role (e.g. airport networks [HI), then the scale- 
free property emerges. When both the two ingredients 
are not neglectable (e.g. bus networks [II]), the degree 
distribution is intervenient between power-law and expo- 
nential ones. 

Fig. 3 shows the simulation results for a = 1 and 
a = 0.5. The degree distribution follows a power-law 
form when a = 1, which well agrees with the analytic 
solution. In the case of a = 0.5, the degree distri- 
bution is more exponential. However, it is remarkably 
broader than that of the Erdos-Renyi model [25[ . Note 
that, the positions of all the nodes are not completely 
uniformly distributed, which will affect the degree distri- 
bution. This effect becomes more prominent when the 
geographical ingredient plays a more important role (i.e. 
smaller a). Therefore, although the simulation result for 
a = 0.5 is in accordance with the analysis qualitatively, 
the quantitative deviation can be clearly observed. 

IV. TOPOLOGICAL DISTANCE 

Denote by the topological distance between the 
ith node and the first node. By using mathematical 
induction, we can prove that there exists a positive 
constant M, such that l(t) < Mln(t). This proposi- 
tion can be easily transferred to prove the inequality 




FIG. 3: (Color online) The degree distributions for the cases 
of a = 1 (a) and a = 0.5 (b). The black squares and red 
curves represent simulated and analytic results, respectively. 
All the data are averaged over 100 independent runs, with 
network size N = 10000 (i.e. t = 10000) fixed. 



l(t + l) < Mln(t + 1) under the condition < Mln(i) 
for 1 < i < t. Indeed, since the network has a tree struc- 
ture, does not depend on time t. Under the frame- 
work of the mean field theory, the iteration equation for 
l(t) reads 

l(t+l) = ^i<»<* ' v > w +lj ( 7 ) 

with the initial condition ^(1) = 0. Eq. (7) can be under- 
stood as follows: At the (t+ l)th time step, the (t+ l)th 
node has probability kf(t)/At to connect with the ith 
node. Since the average topological distance between the 
ith node and the first node is the topological dis- 
tance of the (t + l)th node to the first one is + 1 
if it is connected with the ith node. According to the 
induction assumption, 

l(t+l) < ~ -\ +1- (8) 



4 



Note that, statistically, fej(i) > kj(t) if i < j, therefore 

E K(t)Mi) < E < fcQ W> ln (*) = £ ^ ln «> 

l<i<t \<i<t l<i<t 

(9) 

where (•) denotes the average over all the nodes. Substi- 
tuting inequality (9) into (8), we have 

^+!)<T £ ln « + !- (10) 

l<i<t 

Rewriting the sum in continuous form, we obtain 

M f* 

l(t + l)<— ln(i)da; + 1< Mln(i + 1). (11) 
t Ji 

According to the mathematical induction principle, we 
have proved that the topological distance between the 
«th node and the first node, denoted by could not 
exceed the order 0(ln(i)). For arbitrary nodes i and j, 
clearly, the topological distance between them could not 
exceed the sum + thus the average topological 
distance (d) of the whole network could not exceed the 
order 0(ln(z)) either. This topological characteristic is 
referred to as the small-world effect in network science 
[l|, and has been observed in a majority of real networks. 
Actually, one is able to prove that the order of l(t) in the 
large t limit is equal to ln(i) (see Appendix A for details). 
Furthermore, the iteration equation 

K*+i)= % / y +«(*), as) 

Ei=i /(*) 

for general functions f(i) and a(i), has the following so- 
lution: 

l(t + !)=/(!) + a(t) + (13) 

j=l l^i=l JW 

For the two special cases of a = (a(i) — 1, f(i) = 1) 
and a = 1 (a(i) = 1, /(«) = l/yi), the solutions are 
simply i(t) = ln(t) and l(t) = |ln(t), respectively. 

In Fig. 4, we report the simulation results about 
the average distance (d) vs network size t. In each 
case, the data points can be well fitted by a straight 
line in the semi-log plot, indicating the growth tendency 
(d) ~ ln(t), which agrees well with the analytical solu- 
tion. 



V. EDGE LENGTH DISTRIBUTION 

Denote by ejj the edge between nodes i and j, and 
the geographical length of edge ey is ry = \r% — ¥j\. 
When the (N + l)th node is added to the network, the 
geographical length of its attached edge approximately 
obeys the distribution 

N 2 

Q(r) = 27rr^V 7r % E3fc? , (14) 

i=l 



t ' r 
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FIG. 4: (Color online) The average distance vs network size 
for the cases of a = 1 (a) and a = 0.5 (b). As shown in 
each inset, the data points can be well fitted by a straight 
line in the semi-log plot, indicating the growth of the average 
distance approximately obeys the form (d) ~ ln(t). All the 
data are averaged over 5 independent runs, where the maxi- 
mal network size is N = 10000 (i.e. t = 10000). 



where r <C 1 in the large N limit. The derivation of 
this formula is described as follows. The probability of 
the edge length being between r and r + dr is given by 
the summation ^\ Qidr, where Qidr is the probability 
that ft jv+1 f aus between r and r + dr, and the node i 
minimizes the quantity |rjv+i — ^l 2 /^? among all the 
N previously existing nodes. This probability is approx- 
imately given by 

Qi « ^rH^l-^j 

(15) 
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Straightforwardly, the geographical length distribution of 
the newly added edge at the tth time step (the ith edge 
for short) is obtained as 



P(t,r) = 2vrr^\ 



2irrt \ due-* Atr ' Fa ^ 
Ji/t 

2nrt t due- vAtr2 / Fa ^ u \ 
Jo 



(16) 



The lower boundary in the integral is replaced by in the 
last step, which is valid only when a < 1. The cumulative 
length distribution of the edges at time step T is given 
by 



P(r) 



1 

T 



P(t, r)dt 



irr 3 A 2 T J 



duF 



2a 



-irAr 2 / F° 



nAr 2 

pa 



-irATr /F° 



nATr 2 

pa 



(17) 



where the argument of function F is 1/u. For 1/y/T <C 
r«l, the approximate formula for P{r) reads 



P(r) 



r 3 A 2 T 



duF 2 



and, when r -C 1/VT, 



P(r) rs irrT. 



(18) 



(19) 



If a = 1, the last step in Eq. (16) is invalid but the 
analytic form for P(t, r) can be directly obtained as 



f 1 2 
P{t, r) = 2-nrt / due- 2lxtr *^ 

Jilt 



= — 1=-[(1 + 27rr 2 v^)e- 27rr2 ^ - (1 + 2Trr 2 t)e- 27vr2t ]{20) 
7r r 6 t 

Therefore, when 1/T 1 ^ 4 « r « 1, P(f) is approximately 
given by 



P(r) 



1 .ln. C 



r*T 2nr 2 ' 



(21) 



where C is a numerical constant, and when r <C 1/VT, 
P{r) has the same form as that of Eq. (19). 

Fig. 5 plots the cumulative edge-length distributions. 
From this figure, one can see a good agreement between 
the theoretical and the numerical results. Furthermore, 
one can calculate the expected value of the tth edge's 
geographical length, as 

= f*rP(t,r)dr = _J_ g duF^I 2 {l/u) 
J^P(t,r)dr 2VAi ^duF a {l/u) ' 
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FIG. 5: (Color online) The edge length distributions for the 
cases of a = 1 (a) and a — 0.5 (b). The black squares and red 
curves represent simulation and analytic results, respectively. 
All the data are averaged over 100 independent runs, with 
network size N = 10000 (i.e. T = 10000) fixed. 



which is valid only for sufficiently large t and a < 1. Ac- 
cording to Eq. (22), f(t) decreases as as t increases, 
which is consistent with the intuition since all the t nodes 
are embedded into a 2-dimensional Euclidean space. It 
may also be interesting to calculate the total length R(t) 
of all the edges at the time step t, as 



t 



dt'r{t') 



t ^duF 3a l 2 {\/u) 

A fiduF<*(l/u) ■ 
(23) 

R(t) is proportional to \/t for 1 > a > 0. When a > 1, 
a finite fraction of nodes will be connected with a single 
hub node and therefore we expect that R(t) ~ t in this 
case. Therefore, in the large t limit, R(t) will increase 
quite abruptly when the parameter a exceeds 1. This 
tendency is indeed observed in our numerical simulations, 
as shown in Fig. 6. 
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VI. GEOGRAPHICAL DISTANCE 

For an arbitrary path iq — > i\ — > • • • — > i n from 
node io to i„, the corresponding geographical length 



V™" 1 r- 



where denotes the length of edge 



is 

6y . Accordingly, the geographical distance between two 
nodes is defined as the minimal geographical length of 
all the paths connecting them. Now, we calculate the 
geographical distance s(i) between the ith node and the 
first node. Since our network is a tree graph, s(i) does 
not depend on time. By using the mean field theory, we 
have 



s(t + 1) = — — + r(t + 1), 



or 



s (* + 1) = -7 / d uF a (-) s(ut) + - 7 L. 

where, according to Eq. (22), 

1 [,!: duF 3a / 2 (l/<u) 
W = 



2v A J*duF a (l/u) 



(24) 



(25) 



(26) 



It is not difficult to see that s(t) has an upper bound 
as t approaches infinity. One can use the trial solution 
s(t) = B — C/t" to test this conclusion: 

B - C/(t + if = B - CE/{At f3 ) + W/Vt + l, (27) 
where 



E 



i/' 



duF a {l/u){l/u fs ). 



(28) 



From Eq. (27), one obtains that = 1/2. 

Similar to the solution of Eq. (13), s(t) — > 1 as t — » oo 
for a = Oands(i) -> l/v^ as i — > oo for a = 1. However, 
it only reveals some qualitative property, and the exact 
numbers are not meaningful. This is because the value of 
s(t) is obtained by the average over infinite configurations 
for infinite t, while in one evolving process s(t) is mainly 
determined by the randomly assigned coordinates of the 
tth node. 



VII. CONCLUSION AND DISCUSSION 

In many real-life transportation networks, the geo- 
graphical effect can not be ignored. Some scientists pro- 
posed certain global optimal algorithms to account for 
the geographical effect on the structure of a static net- 
work [19L |2(J ] . On the other hand, many real networks 
grow continuously. Therefore, we proposed a growing 
network model based on an optimal policy involving both 
topological and geographical measures. We found that 
the degree distribution will be broader when the topolog- 
ical ingredient plays a more important role (i.e. larger a), 




FIG. 6: The total edge length R(t) vs parameter a with 
t = 10000 fixed. R(t) sharply increases when a exceeds the 
critical value a c = 1, which well agrees with the theoretical 
prediction. 



and when a exceeds a critical value a c — 1, a finite frac- 
tion of nodes will be connected with a single hub node and 
the geographical effect will become insignificant. This 
critical point can also be observed when detecting the 
total geographical edge-length R(t) in the large t limit. 
We obtained some analytical solutions for degree distri- 
bution, edge-length distribution, and topological as well 
as geographical distances, based on reasonable approxi- 
mations, which are well verified by simulations. 

Although the present model is based on some ideal as- 
sumptions, it can, at least qualitatively, reproduce some 
key properties of the real transportation networks. In 
Table 1, we list some empirical degree distributions of 
transportation networks. Clearly, when building a new 
airport, we tend to firstly open some flights connected 
with previously central airports which are often of very 
large degrees. Even though the central airports may be 
far from the new one, to open a direct flight is relatively 
convenient since one doesn't need to build a physical link. 
Therefore, the geographical effect is very small in the ar- 
chitecture of airport networks, which corresponds to the 
case of larger a that leads to an approximately power- 
law degree distribution. For other four cases shown in 
Table 1, a physical link, which costs much, is necessary 
if one wants to connect two nodes, thus the geographical 
effect plays a more important role, which corresponds to 
the case of smaller a that leads to a relatively narrow 
distribution. 

A specific measure of geographical network is its edge- 
length distribution. A very recent empirical study [191 ] 
shows that the edge- length distribution of the highly het- 
erogenous networks (e.g. airport networks, correspond- 
ing to the present model with larger a) displays a single- 
peak function with the maximal edge-length about five 
times longer than the peaked value (see Fig. lc of Ref. 
[l||), while in the extreme homogenous networks (e.g. 
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TABLE I: Empirical degree distributions for geographical 
transportation networks. 



Data-Class 



Data-Set 



Degree Distribution 



Power Grid 



Southern California [26] 
Whole US [15] 
Italian [141 



Exponential™ 
p(k) 
P(k) ~ 



-0.55fc 



Subway 



Boston [28] 
Seoul [29] 
Tokyo [29] 



Narrow: fc max < 10 
Narrow: fc max < 10 
Narrow: fc max < 10 



Railway 



Indian [10] 
Switzerland [30] 
Central Europe [301 



p(k) ~ e 

Exponential 
Exponential 



Bus & Tramway 



Krakow [311 
Warsaw [31] 
Szczecin [31] 
Bialystok [31] 



Narrow: fc max = 11 
Narrow: fc max — 13 
Narrow: fc max = 18 
Narrow: fc max = 19 



Airport 



World-Wide [11] 
Chinese [32] 
Indian [33] 



p(fc) ~ k~ 2 - Uc 
p(fe) ~fc- 2 ' 05d 
p(Jfe) ~ k- 2/2 



"In Ref. 0, the authors claimed that this distribution follows a 
power-law form p(k) ~ k ' with exponent 7 4. Actually, when 
7 gets larger, the network will topologically gets closer to random 
graph [23] - 

''In Ref. [30ll . the authors displayed an exponential degree distri- 
bution of mass transportation networks in Warsaw. 
c Actually a truncated power-law distribution. 
"Actually a double power-low distribution. 



railway networks, corresponding to the present model 
with a — > 0), only the very short edge can exist (see Fig. 
la of Ref. [Ijl]). These empirical results agree well with 
the theoretical predictions of the present model. Firstly 
when a is obviously larger than zero, the edge-length dis- 
tribution is single-peaked with its maximal edge-length 
about six times longer than the peak value (see Fig. 5). 
And, when a is close to zero, Eq. (17) degenerates to the 
form 



P(r) = 



(29) 



where N denotes the network size. Clearly, in the large 
TV limit, except a very few initially generated edges, only 
the edge of very small length r can exist. 

The analytical approach is only valid for the tree struc- 
ture with m = 1. However, we have checked that all the 
results will not change qualitatively if m is not too large 
compared with the network size. Some analytical meth- 
ods proposed here are simple but useful, and may be ap- 
plied to some other related problems about the statistical 
properties of complex networks. For example, a similar 
(but much simpler) approach, taken in section 4, can also 
be used to estimate the average topological distance for 
some other geographical networks [3J, [35| . 

Finally, it is worthwhile to emphasize that, the geo- 
graphical effects should also be taken into account when 
investigating the efficiency (e.g. the traffic throughput 
(36|) of transportation networks. Very recently some au- 
thors started to consider the geographical effects on dy- 
namical processes, such as epidemic spreading 37J and 



cascading . over scale-free networks. We hope the 
present work can further enlighten the readers on this 
interesting subject. 
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APPENDIX A: THE SOLUTION OF l(t) 



Substituting Eq. (4) into Eq. (7), one obtains that 



«(* + !) = 4/ duF a (- ) l(ut) + 1. 



Then, define 



B Ail duF a (l/u)\n(l/u) 



(Al) 



(A2) 



We next prove that l(t) w Bint in the large t limit by 
using mathematical induction. Suppose for sufficiently 
large t, all are less than Cln(z) for i < t with C 
being a constant greater than B. Then, from Eq. (Al), 
we have 



Ut+l) < ^ f duF a ( - 
~ Aj 1/t \u 



ln(ut) + 1 



= ~ + l + M\n{t)+o(j 
< Cln(t + 1) 



(A3) 



Therefore, < Cln(«) for all i. Similarly, suppose 
for sufficiently large t, all are greater than Qhii for 
i < f with Q being a constant less than B. Then, from 
Eq. (Al), we have 



~ AJ 1/t \u 



ln(ut) + 1 



= -| + l + Mln(t)+o(j 
> Qln(t + 1). 



(A4) 



Therefore, > Qln(i) for all i. 

Combine both the upper bound (A3) and lower bound 
(A4), we obtain the order of l(t) in the large t limit, as 
l(t) «51n(t). 
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